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ABSTRACT 

Using  the  Navy  Precision  Optical  Interferometer,  we  measured  the  angular  diameters  of  10  stars  that  have  previously 
measured  solar-like  oscillations.  Our  sample  covered  a  range  of  evolutionary  stages  but  focused  on  evolved  subgiant 
and  giant  stars.  We  combined  our  angular  diameters  with  Hipparcos  parallaxes  to  determine  the  stars’  physical  radii, 
and  used  photometry  from  the  literature  to  calculate  their  bolometric  fluxes,  luminosities,  and  effective  temperatures. 

We  then  used  our  results  to  test  the  scaling  relations  used  by  asteroseismology  groups  to  calculate  radii  and  found 
good  agreement  between  the  radii  measured  here  and  the  radii  predicted  by  stellar  oscillation  studies.  The  precision 
of  the  relations  is  not  as  well  constrained  for  giant  stars  as  it  is  for  less  evolved  stars. 

Key  words:  stars:  fundamental  parameters  -  stars:  oscillations  -  techniques:  high  angular  resolution  -  techniques: 
interferometric 


1.  INTRODUCTION 

Asteroseismology,  the  study  of  stellar  oscillations,  is  a 
powerful  tool  to  infer  information  about  stellar  structure  with 
minimal  model  dependence  (see,  e.g..  Brown  &  Gilliland  1994; 
Christensen-Dalsgaard  2004).  The  frequencies  of  the  observed 
oscillations  depend  on  the  sound  speed  inside  the  star,  which  in 
turn  is  dependent  on  properties  of  the  interior  such  as  density, 
temperature,  and  gas  motion  (Carrier  et  al.  2010).  The  num¬ 
ber  of  stars  observed  using  asteroseismology  and  the  quality 
of  the  data  have  increased  dramatically  in  recent  years,  due  to 
the  photometric  space  missions  MOST  ( Microvariability  and 
Oscillations  ofSTars;  Walker  et  al.  2003),  CoRoT  ( Convection , 
Rotation,  and  planetary  Transits',  Baglin  et  al.  2006;  Auvergne 
et  al.  2009),  and  Kepler  (Borucki  et  al.  2010;  Koch  et  al.  2010). 
The  resulting  stellar  parameters  are  key  to  the  statistical  analysis 
of  fundamental  stellar  properties  and  for  testing  stellar  interior 
and  evolutionary  models  (see,  e.g.,  Chaplin  et  al.  201 1). 

Interferometry  has  the  potential  to  make  important  contribu¬ 
tions  to  asteroseismology,  in  part  through  the  determinations 
of  the  targets’  sizes  (Cunha  et  al.  2007).  Using  interferometry, 
we  can  measure  the  angular  diameters  of  stars  with  resolutions 
down  to  tenths  of  a  milliarcsecond  (see,  e.g.,  Huber  et  al.  2012a; 
Baines  et  al.  2012).  Once  we  know  the  apparent  diameter  of 
a  star  as  well  as  its  distance  from  parallax  measurements,  we 
can  calculate  its  physical  size.  Then  we  can  test  the  relation¬ 
ships  used  to  derive  stellar  properties  from  asteroseismology 
observations  by  comparing  the  radii  estimated  using  the  astero¬ 
seismology  relations  to  those  measured  interferometrically. 

Huber  et  al.  (2012b)  presented  interferometric  diameters  of 
10  stars  that  had  oscillation  measurements  from  CoRoT  and 
Kepler.  They  found  an  agreement  between  asteroseismic  and 
interferometric  radii  of  <4%  for  dwarf  stars  and  ~13%  for 
giant  stars.  Their  sample  included  five  dwarf,  one  subgiant,  and 
four  giant  stars.  Here  we  focus  on  the  more  evolved  stars:  one 
dwarf,  four  subgiant,  and  five  giant  stars. 

We  observed  these  stars  using  the  Navy  Precision  Optical  In¬ 
terferometer  (NPOI)  in  order  to  measure  their  angular  diameters. 
We  then  calculated  their  radii  and  effective  temperatures,  and 
used  spectral  energy  distribution  fits  to  determine  their  bolomet¬ 
ric  fluxes  and  luminosities.  Section  2  discusses  the  NPOI  and 


our  observing  process;  Section  3  describes  the  visibility  mea¬ 
surements  and  how  we  calculated  various  stellar  parameters; 
Section  4  explores  the  relationship  between  radii  determined 
using  asteroseismology  observations  and  radii  measured  inter¬ 
ferometrically;  and  Section  5  summarizes  our  findings. 

2.  INTERFEROMETRIC  OBSERVATIONS 

The  NPOI  is  an  interferometer  located  on  Anderson  Mesa, 
AZ,  and  consists  of  two  nested  arrays:  the  four  stations  of 
the  astrometric  array  (AC,  AE,  AW,  and  AN,  which  stand 
for  astrometric  center,  east,  west,  and  north,  respectively)  and 
the  six  stations  of  the  imaging  array,  of  which  two  stations 
are  currently  in  operation  (E6  and  W7)  and  three  more  will 
be  coming  online  in  the  near  future  (E7,  E10,  and  W10).  The 
current  baselines,  i.e.,  the  distances  between  the  stations,  range 
from  1 6  to  79  m,  and  our  maximum  baseline  will  be  432  m  when 
the  E10  and  W10  stations  are  completed  within  the  next  year. 
We  use  a  12  cm  region  of  the  50  cm  siderostats  and  observe 
in  16  spectral  channels  spanning  550-850  nm  simultaneously 
(Armstrong  et  al.  1998). 

Each  observation  consisted  of  a  30  s  coherent  (on  the 
fringe)  scan  in  which  the  fringe  contrast  was  measured  every 
2  ms,  paired  with  an  incoherent  (off  the  fringe)  scan  used  to 
estimate  the  additive  bias  affecting  the  visibility  measurements 
(Hummel  et  al.  2003).  Scans  were  taken  on  up  to  three  baselines 
simultaneously.  Each  coherent  scan  was  averaged  to  1  s  data 
points,  and  then  to  a  single  30  s  average.  The  dispersion  of  1  s 
points  provided  an  estimate  of  the  internal  uncertainties. 

The  target  list  was  derived  from  the  sample  of  stars  with 
stellar  oscillations  that  were  bright  enough  to  observing  using 
the  NPOI,  which  has  a  magnitude  limit  of  V  =  6.5.  They  also 
had  to  be  resolved  with  the  longest  existing  baseline,  which 
gives  a  resolution  limit  of  approximately  1  mas.  This  resulted  in 
a  list  of  10  targets  with  stellar  oscillation  observations  available 
to  observe  using  the  NPOI. 

We  interleaved  data  scans  of  the  10  asteroseismic  targets  with 
one  to  three  calibrator  stars  for  each  target.  Our  calibrators  are 
stars  that  are  significantly  less  resolved  on  the  baselines  used 
than  the  targets.  This  meant  that  uncertainties  in  the  calibrator’s 
diameter  did  not  affect  the  target’s  diameter  calculation  as  much 
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as  if  the  calibrator  star  had  a  substantial  angular  size  on  the  sky. 
The  calibrator  and  target  scans  were  measured  as  close  in  time 
and  space  as  possible,  which  allowed  us  to  convert  instrumental 
target  and  calibrator  visibilities  to  calibrated  visibilities  for 
the  target.  Preference  was  given  to  calibrators  within  10°  of 
the  target  stars,  as  was  the  case  for  13  of  the  16  calibrator  stars 
used.  On  rare  occasions,  no  suitable  calibrator  stars  were  within 
that  angular  distance  so  we  resorted  to  stars  that  were  more 
distant,  with  a  maximum  separation  of  17°. 

We  estimated  the  calibrator  stars’  sizes  by  constructing  their 
spectral  energy  distribution  (SED)  fits  using  photometric  val¬ 
ues  published  in  Ljunggren  &  Oja  (1965),  McClure  &  Forrester 
(1981),  Olsen  (1993),  Jasevicius  et  al.  (1990),  Golay  (1972), 
Haggkvist  &  Oja  (1970),  Kornilov  et  al.  (1991),  Eggen  (1968), 
Johnson  et  al.  (1966),  Cutri  et  al.  (2003),  and  Gezari  et  al. 
(1993)  as  well  as  spectrophotometry  from  Glushneva  et  al. 
(1983,  1998),  and  Kharitonov  et  al.  (1997)  obtained  via  the 
interface  created  by  Mermilliod  et  al.  (1997).  The  assigned  un¬ 
certainties  for  the  Two  Micron  All  Sky  Survey  infrared  measure¬ 
ments  are  as  reported  in  Cutri  et  al.  (2003),  and  an  uncertainty 
of  0.05  mag  was  assigned  to  the  optical  measurements.  We  de¬ 
termined  the  best  fit  stellar  spectral  template  to  the  photometry 
from  the  flux-calibrated  stellar  spectral  atlas  of  Pickles  (1998) 
using  the  /2  minimization  technique  (Press  et  al.  1992;  Wall  & 
Jenkins  2003).  The  resulting  calibrator  angular  diameter  esti¬ 
mates  are  listed  in  Table  1 . 


3.  RESULTS 


3.1.  Angular  Diameter  Measurement 

Interferometric  diameter  measurements  use  V2,  the  square  of 
the  fringe  visibility.  For  a  point  source,  V2  is  unity,  while  for 
a  uniformly-illuminated  disk,  V2  —  [2J\(x)/x]2,  where  J\  is 
the  Bessel  function  of  the  first  order,  x  =  n  BO^n'k  1 ,  B  is  the 
projected  baseline  toward  the  star’s  position,  6*ud  is  the  apparent 
uniform  disk  angular  diameter  of  the  star,  and  k  is  the  effective 
wavelength  of  the  observation  (Shao  &  Colavita  1992).  0\jd 
results  are  listed  in  Table  2.  Our  data  files  in  OIFITS  format  are 
available  upon  request. 

A  more  realistic  model  of  a  star’s  disk  includes  limb  darkening 
(LD).  If  a  linear  LD  coefficient  px  is  used, 


V2 


(1  -  ah) 


A(xld) 

xld 


1/2  73/2 Old) 


X 


3/2 

LD 


(1) 


where  .vld  =  tt  B0\  n'D  1  (Hanbury  Brown  et  al.  1974).  We  used 
effective  temperature  (Teff)  and  surface  gravity  (log  g)  values 
from  the  literature  with  a  microturbulent  velocity  of  2  km  s'1 
and  to  obtain  /x;  from  Claret  &  Bloemen  (2011).  These  values 
and  the  resulting  Old  are  listed  in  Table  2.  Figures  1  and  2 
show  the  9ld  fits  for  all  the  stars.  The  two  stars  with  the  largest 
percent  uncertainties  in  the  Old  fit  (2%)  are  HD  146791  and 
HD  181907.  This  is  because  their  visibility  curves  are  less  well 
sampled  with  respect  to  spatial  frequency  than  the  other  targets. 
The  remaining  stars  have  uncertainties  in  Old  of  0.2%  to  0.6%. 

Of  the  10  stars  measured  here,  7  had  previous  interferometric 
diameter  measurements.  They  are  listed  in  Table  3  and  are 
plotted  against  our  values  in  Figure  3.  In  all  cases  but  one, 
the  uncertainty  on  our  diameter  measurement  is  smaller  than 
those  from  the  literature,  and  they  all  agree  to  within  3 a . 


The  uncertainty  for  the  Old  fit  was  derived  using  the  method 
described  in  Tycner  et  al.  (2010),  who  showed  that  a  non-linear 
least-squares  method  does  not  sufficiently  account  for  atmo¬ 
spheric  effects  on  time  scales  shorter  than  the  window  between 
target  and  calibrator  observations.  They  instead  describe  a  boot¬ 
strap  Monte  Carlo  method  that  treats  the  observations  as  groups 
of  data  points  because  the  NPOI  collects  data  in  scans  con¬ 
sisting  of  16  channels  simultaneously.4  They  discovered  that 
when  the  16  data  points  were  analyzed  individually,  a  single 
scan’s  deviation  from  the  trend  had  a  large  impact  on  the  re¬ 
sulting  diameter  and  uncertainty  calculation.  On  the  other  hand, 
when  they  preserved  the  inherent  structure  of  the  observational 
data  by  using  the  groups  of  16  channels  instead  of  individual 
data  points,  the  uncertainty  on  the  angular  diameter  was  larger 
and  more  realistic.  This  method  makes  no  assumptions  about 
underlying  uncertainties  due  to  atmospheric  effects,  which  are 
applicable  to  all  stars  observed  using  ground-based  instruments. 
It  should  be  noted  that  the  number  of  calibrator  stars  used  in  the 
observations  have  no  apparent  effect  on  the  Old  fit  uncertainty. 

3.2.  Stellar  Radius,  Luminosity,  and  Effective  Temperature 

For  each  star,  the  parallax  from  van  Leeuwen  (2007)  was 
converted  into  a  distance,  which  we  then  combined  with  our 
measured  Old  to  calculate  the  linear  radius  ( R ).  In  order  to 
determine  the  luminosity  ( L )  and  Te .g,  we  constructed  each  star’s 
SED  using  the  sources  and  technique  of  fitting  spectral  templates 
to  observed  photometry  as  described  in  Section  2.  The  resulting 
SED  gave  us  the  bolometric  flux  (Frol)  and  allowed  for  the 
calculation  of  extinction  Ay  with  the  wavelength-dependent 
reddening  relations  of  Cardelli  et  al.  (1989). 

We  combined  our  Fbol  values  with  the  stars’  distances  to 
estimate  L  using  L  =  And1  /’hoi.-  We  also  combined  the  Fbol 
with  Old  to  determine  each  star’s  effective  temperature  by 
inverting  the  relation, 


^BOL  =\elDoT*e,  (2) 

where  a  is  the  Stefan-Boltzmann  constant  and  Old  is  in  radians. 

Because  /x;  is  chosen  based  on  a  given  7’cn  ,  we  checked  to  see 
if  p,x  and  therefore  Old  would  change  based  on  our  new  Teg.  In 
most  cases,  /x>  changed  by  0.0  or  0.01,  and  the  largest  difference 
was  0.08  for  HD  181907.  The  resulting  Old  values  changed  at 
most  by  1%,  and  all  but  three  changed  by  0.2%  or  less.  This  was 
well  within  the  uncertainties  on  Old,  and  re-calculating  7jn  with 
the  new  0LD  made  at  most  a  26  K  difference  (for  HD  181907, 
which  has  an  uncertainty  of  199  K).  These  values  all  converged 
after  this  one  iteration,  and  these  are  the  final  numbers  listed  in 
Table  2. 


4.  DISCUSSION 

Two  scaling  equations  relate  observed  asteroseismic  quanti¬ 
ties  to  fundamental  stellar  parameters: 

An  oc  MiR~T  (3) 

where  Av  is  the  large  separation  of  oscillation  modes  of  the 
same  degree  and  consecutive  orders  and  M  is  the  mass  of  the 


4  For  every  “scan,”  30  s  of  data  are  collected  in  each  of  the  16  wavelength 
channels  with  a  measurement  once  every  2  ms.  During  the  processing 
described  in  Section  2,  all  the  30  s  scan’s  data  points  are  averaged  into  one  data 
point  for  each  channel,  so  we  go  from  30  s  of  data  per  channel  to  one  averaged 
data  point  per  channel. 
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Table  1 

Observing  Log  and  Calibrator  Stars’  Angular  Diameters 


Target 

HD 

Other 

Name 

Calibrator 

HD 

Date 

(UT) 

Baselines 

Used3 

No. 

Obs 

#LD,cal 

(mas) 

10700 

r  Cet 

11171 

2005  Sep  17 

AC-W7,  AE-W7 

19 

0.682  ±  0.034 

2005  Sep  18 

AC-W7 

9 

2005  Sep  25 

AC-W7,  AE-W7 

44 

121370 

1 1  Boo 

122408 

2012  May  10 

AW-E6,  E6-W7 

126 

0.521  ±  0.026 

2012  May  11 

AE-AW,  E6-W7 

90 

2012  May  12 

E6-W7 

27 

2012  May  15 

AW-E6,  E6-W7 

108 

146791 

6  Oph 

141513 

2010  Mar  21 

AW-W7 

30 

0.524  ±  0.026 

2010  Mar  28 

AW-W7 

13 

150680 

£  Her 

156164 

2005  Aug  18 

AE-W7,  E6-W7 

50 

0.887  ±  0.044 

2005  Aug  25 

AC-W7,  AE-W7,  E6-W7 

120 

153210 

k  Oph 

148112 

2013  Feb  23 

AC-AE,  AC-AW,  AC-E6 

18 

0.443  ±  0.022 

2013  Mar  1 

AC-AW 

12 

2013  Mar  2 

AC-E6,  AW-E6 

63 

2013  Mar  5 

AC-E6,  AW-E6 

198 

2013  Mar  6 

AC-E6,  AW-E6 

72 

2013  Mar  7 

AC-E6,  AW-E6 

90 

152614 

2013  Feb  23 

AC-E6,  AC-AW 

34 

0.316  ±0.016 

2013  Feb  25 

AC-E6 

14 

2013  Mar  2 

AC-E6,  AW-E6 

90 

2013  Mar  5 

AC-E6,  AW-E6 

261 

2013  Mar  7 

AC-E6,  AW-E6 

63 

2013  Apr  1 

AE-E6,  AW-E6 

53 

2013  Apr  2 

AE-E6,  AW-E6 

35 

2013  Apr  3 

AE-E6,  AW-E6 

54 

147547 

2013  Apr  19 

AE-AW,  AW-E6 

68 

0.970  ±  0.049 

2013  Apr  22 

AE-E6,  AW-E6 

115 

2013  Apr  23 

AE-E6,  AW-E6 

91 

2013  Apr  26 

AE-E6,  AW-E6 

165 

161797 

/i  Her 

166014 

2010  Aug  21 

AE-AW,  AW-E6 

270 

0.596  ±  0.030 

2010  Aug  26 

AE-AW,  AW-E6 

58 

163588 

|  Dra 

159541 

2013  Apr  27 

AE-AW,  AW-E6 

18 

0.520  ±  0.026 

2013  Apr  28 

AW-E6 

43 

2013  Apr  29 

AE-AW,  AW-E6 

97 

168151 

2013  Apr  27 

AE-AW,  AW-E6 

43 

0.655  ±  0.033 

2013  Apr  28 

AE-AW,  AW-E6 

49 

2013  Apr  29 

AE-AW,  AW-E6 

90 

2013  May  1 

AW-E6 

10 

2013  May  3 

AE-AW 

12 

2013  May  9 

AW-E6 

170 

2013  May  12 

AW-E6 

23 

2013  May  13 

AW-E6 

236 

2013  May  14 

AE-AW 

32 

184006 

2013  May  1 

AW-E6 

32 

0.703  ±  0.035 

2013  May  12 

AE-AW,  AW-E6 

48 

2013  May  13 

AW-E6 

104 

2013  May  14 

AE-AW 

49 

168723 

rj  Ser 

161868 

2004  May  20 

AC-AE,  AC-AW 

90 

0.668  ±  0.033 

2004  May  24 

AC-AE,  AC-AW 

17 

164353 

2007  May  18 

AW-E6,  E6-W7 

120 

0.371  ±  0.022 

2007  May  24 

AN-E6,  AW-E6,  E6-W7 

278 

2007  May  25 

AN-E6,  AW-E6 

20 

2007  May  31 

AW-E6,  E6-W7 

128 

181907 

HR  7349 

177756 

2013  May  21 

AW-E6 

62 

0.571  ±  0.029 

2013  May  31 

AC-AW,  AW-E6 

16 

2013  Jun  3 

AW-E6 

23 

2013  Jun  5 

AW-E6 

10 

184930 

2013  May  21 

AW-E6 

71 

0.311  ±0.016 

2013  May  31 

AC-AW,  AW-E6 

16 

2013  Jun  3 

AW-E6 

21 

2013  Jun  5 

AW-E6 

10 

188512 

P  Aql 

195810 

2007  May  26 

AN-AW,  AW-W7 

14 

0.394  ±  0.020 

2007  May  30 

AW-E6,  E6-W7 

57 

2007  May  31 

AN-E6,  AW-E6,  E6-W7 

219 

2007  Jun  4 

AW-E6,  E6-W7 

140 

2007  Jun  9 

AN-E6,  AW-E6,  E6-W7 

120 

2007  Jun  11 

AW-E6,  E6-W7 

89 

Notes.  a  The  maximum  baseline  lengths  are  AC-AE  18.9  m,  AC- AW  22.2  m,  AC-E6  34.4  m,  AC-W7  51.3  m,  AE-AW  37.5  m,  AE-W7  64.2  m,  AN-AW  38.2  m,  AN-E6  45.6  m, 
AW-E6  53.3  m,  AW-W7  29.5  m,  and  E6-W7  79.4  m.  The  0LD,cai  estimates  were  determined  using  the  technique  described  in  Section  2. 
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Figure  1.  Old  fits  for  stars  observed  with  one  calibrator.  The  solid  lines  represent  the  theoretical  visibility  curve  for  the  best  fit  Oli ),  the  points  are  the  calibrated 
visibilities,  and  the  vertical  lines  are  the  measurement  uncertainties.  The  uncertainty  in  the  Old  fit  is  not  shown  because  it  largely  indistinguishable  from  the  best  fit 
Old  curve  on  this  scale.  See  Table  2  for  the  uncertainty. 


Table  2 

Stellar  Parameters 


Target 

HD 

Spectral 

Type 

Parallax 

(mas) 

Mx 

#UD 

(mas) 

Old 

(mas) 

&LD 

(%) 

No. 

Cals 

R 

(Rq) 

L 

(Le) 

Fbol 

(10~8  erg  s-1  cm^2) 

Tts 

(K) 

tTTeff 

% 

10700 

G8.5  V 

273.96  ±0.17 

0.60 

1.952  ±0.003 

2.072  ±  0.010 

0.5 

i 

0.81  ±0.01 

0.47  ±  0.01 

113.0  ±0.3 

5301  ±  13 

0.2 

121370 

GO  IV 

87.75  ±  1.24 

0.52 

2.023  ±  0.002 

2.134  ±  0.012 

0.6 

i 

2.61  ±0.04 

8.69  ±  0.25 

214.0  ±0.7 

6128  ±  18 

0.3 

146791 

G9.5  III 

30.64  ±  0.20 

0.65 

2.772  ±  0.007 

2.966  ±  0.061 

2.1 

i 

10.40  ±  0.22 

56.64  ±  1.29 

170.0  ±  3.2 

4907  ±  55 

1.1 

150680 

GO  IV 

93.32  ±  0.47 

0.53 

2.175  ±0.001 

2.266  ±  0.014 

0.6 

i 

2.61  ±0.02 

8.12  ±0.08 

226.0  ±  0.6 

6029  ±  19 

0.3 

153210 

K2  III 

35.66  ±  0.20 

0.71 

3.479  ±  0.001 

3.657  ±  0.013 

0.4 

3 

11.02  ±0.07 

50.67  ±  1.14 

206.0  ±  4.0 

4367  ±  24 

0.5 

161797 

G5IV 

123.33  ±  0.16 

0.58 

1.851  ±0.002 

1.957  ±0.012 

0.6 

1 

1.71  ±0.01 

2.10  ±0.01 

102.0  ±  0.2 

5317  ±  16 

0.3 

163588 

K2  III 

28.98  ±0.12 

0.70 

2.894  ±  0.002 

3.116  ±0.008 

0.3 

3 

11.56  ±0.06 

47.30  ±  0.44 

127.0  ±  0.6 

4451  ±  7 

0.2 

168723 

K0  III  IV 

53.93  ±0.18 

0.65 

2.852  ±  0.001 

2.970  ±  0.007 

0.2 

2 

5.92  ±  0.02 

17.85  ±0.13 

166.0  ±  0.4 

4875  ±  7 

0.1 

181907 

G8  III 

9.67  ±  0.34 

0.67 

1.038  ±0.009 

1.089  ±0.023 

2.1 

2 

12.10  ±0.50 

96.68  ±  16.0 

29.5  ±  4.3 

5227  ±  199 

3.8 

188512 

G9.5  IV 

78.00  ±  0.20 

0.63 

2.042  ±  0.001 

2.166  ±0.009 

0.4 

1 

2.98  ±0.01 

4.99  ±  0.03 

97.1  ±0.2 

4992  ±  1 1 

0.2 

Notes.  The  parallaxes  are  from  van  Leeuwen  (2007);  the  /ix  coefficients  are  from  Claret  &  Bloemen  (2011)  in  the  J?-band  with  a  microturbulent  velocity  of 
2  km  s-1.  The  sources  of  Teff  and  log  g  used  to  determine  fix  were  the  following:  Prugniel  et  al.  (2011)  for  HD  10700,  HD  121370,  HD  161797,  HD  168723, 
and  HD  188512;  Wu  et  al.  (201 1)  for  HD  146791  and  HD  163588;  Prugniel  et  al.  (2007)  for  HD  150680;  Morel  &  Miglio  (2012)  for  HD  153210;  and  Ammons 
et  al.  (2006)  for  HD  181907. 


star  (Ulrich  1986),  and 

tw  cx  MR-2T;f{\  (4) 

where  vmax  is  the  frequency  of  maximum  oscillation  power 
(Brown  et  al.  1991;  Kjeldsen  &  Bedding  1995).  These  equa¬ 
tions  are  often  used  to  calculate  stellar  radii  and  masses  from 
oscillation  observations.  However,  when  R  is  measured  interfer- 
ometrically,  we  can  test  the  relations  themselves. 

We  used  Ay  and  vmax  from  the  references  listed  in  Table  4 
and  assumed  uncertainties  of  1%  in  Av  and  3%  in  ymax  (Huber 


et  al.  2012b)  when  no  uncertainties  were  provided  in  the 
references.  We  combined  the  frequency  measurements  with 
effective  temperatures  from  the  literature  (7ut)  to  calculate  R 
from  the  asteroseismic  measurements  alone.  It  should  be  noted 
that  the  Tut  used  has  little  impact  on  vmax.  A  variation  of 
100  K  causes  a  0.9%  change  in  vmax,  which  is  typically  on 
the  order  of  or  smaller  than  the  uncertainties  in  7)^  (Huber 
et  al.  2012b). 

Table  4  and  Figure  4  show  the  results  comparing  the  radii 
calculated  from  asteroseismology  (Rd)  and  those  measured  using 
interferometry  (/?;).  The  stars  with  the  largest  difference  between 
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Figure  2.  #ld  fits  for  stars  observed  with  two  or  three  calibrators.  The  symbols  are  the  same  as  in  Figure  1. 


Table  3 

Angular  Diameter  Comparison 


Target 

HD 

$LD,  this work 

(mas) 

$LD,SED 

(mas) 

% 

Diff 

^LD,  previous 

(mas) 

Reference 

% 

Diff 

10700 

2.072  ±  0.010 

2.047  ±  0.038 

1 

1.97  ±0.05 

Pijpers  et  al.  (2003) 

5 

2.015  ±  0.011 

di  Folco  et  al.  (2007) 

3 

121370 

2.134  ±0.012 

2.280  ±  0.069 

7 

2.28  ±  0.07 

Nordgren  et  al.  (2001) 

7 

2.269  ±  0.025 

Mozurkewich  et  al.  (2003) 

6 

2.200  ±  0.027 

Thevenin  et  al.  (2005) 

3 

146791 

2.966  ±  0.061 

3.031  ±0.165 

2 

2.961  ±  0.007 

Mazumdar  et  al.  (2009) 

0.2 

150680 

2.266  ±  0.014 

2.342  ±  0.071 

3 

2.49  ±  0.09 

Nordgren  et  al.  (2001) 

10 

2.367  ±0.051 

Mozurkewich  et  al.  (2003) 

4 

153210 

3.657  ±0.013 

3.960  ±0.199 

8 

N/A 

N/A 

- 

161797 

1.957  ±0.012 

1.767  ±0.051 

10 

1.953  ±  0.039 

Mozurkewich  et  al.  (2003) 

0.2 

1.975  ±  0.025 

Absil  et  al.  (2013) 

1 

163588 

3.116  ±0.008 

3.114  ±0.154 

0.1 

N/A 

N/A 

- 

168723 

2.970  ±  0.007 

2.993  ±0.160 

1 

2.944  ±  0.010 

Merand  et  al.  (2010) 

1 

181907 

1.089  ±0.023 

1.185  ±0.112 

8 

N/A 

N/A 

- 

188512 

2.166  ±0.009 

1.915  ±0.054 

12 

2.18  ±0.09 

Nordgren  et  al.  (1999) 

1 

Notes.  The  SED  fit  was  created  using  the  method  described  in  Section  3.2.  If  more  than  one  diameter  is  available  in  the  literature,  we  used  the  most  recent  one 
when  plotting  the  results  in  Figure  3. 


the  two  are  HD  153210,  HD  161797,  and  HD  168723.  For  the 
latter  two  stars,  previously  published  angular  diameters  agree 
with  our  measurement  to  1%  or  less.  HD  153210  has  not 
been  previously  measured,  and  our  diameter  agrees  with  that 
predicted  by  the  SED  within  2a.  For  the  remaining  stars,  Rd  and 
Ri  agree  within  3a  (~10%)  in  all  cases  except  for  HD  146791 
and  HD  150680,  which  agree  within  4a  (~15%). 

Barban  et  al.  (2004)  quote  a  range  for  vmax  between  80  and 
170  /xHz  for  HD  168723,  and  Stello  et  al.  (2009)  lists  a  value 
of  1 30 //Hz.  If  we  use  130 //Hz  in  our  calculation  for  Rd,  the 


result  is  12.08  ±  0.44  Re,  which  is  approximately  twice  the 
value  of  R[  =  5.92  ±  0.02  Re.  However,  if  we  use  the  lower  end 
of  the  range,  i.e.,  80  /xHz,  Rd  is  7.43  ±  0.27  RQ,  which  is  still  a 
26%  difference  from  the  one  presented  here.  We  believe  further 
asteroseismic  observations  of  this  star  would  be  particularly 
interesting. 

Some  stars  have  more  evenly  and  completely  sampled  data 
along  the  visibility  curves  than  others;  for  example,  HD  146791 
and  HD  181 907  do  not  have  as  wide  a  range  of  measurements 
as  a  function  of  spatial  frequency  that  other  stars  such  as 
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Figure  3.  Comparison  between  #ld  measured  here  and  previous  interferometric  measurements  from  the  literature  (left  panel)  and  compared  to  SED  fits  (right  panel). 
The  bottom  panels  show  the  residuals  to  the  fit.  The  values  used  are  listed  in  Table  3. 


Table  4 

Comparing  Interferometric  and  Asteroseismic  Results 


Target 

HD 

7j 

(K) 

Tit 

(K) 

Ri 

(«o) 

Ra 

(R0) 

% 

diff 

Av 

(MHz) 

t-’max 

(MHz) 

Reference 

M 

( M0 ) 

^max,calc 

(MHz) 

% 

diff 

10700 

5301  ±  13 

5348  ±  45 

0.81  ±  0.01 

0.90  ±  0.03 

12 

169 

4500 

Teixeira  et  al.  (2009) 

0.84  ±  0.02 

4062  ±  108 

10 

121370 

6128  ±  18 

5967  ±  45 

2.61  ±  0.04 

2.86  ±0.10 

9 

39.9  ±0.1 

750 

Carrier  et  al.  (2005) 

1.56  ±0.07 

677  ±  37 

10 

146791 

4907  ±  55 

4918  ±  28 

10.40  ±  0.22 

11.76  ±0.43 

13 

5.3  ±0.1 

60 

Barban  et  al.  (2007) 

1.74  ±0.13 

53  ±  5 

11 

150680 

6029  ±  19 

5758  ±  81 

2.61  ±  0.02 

3.05  ±0.11 

17 

37.01 

701a 

Kallinger  et  al.  (2009) 

1.34  ±0.04 

586  ±21 

16 

153210 

4637  ±  24 

4559  ±  116 

11.02  ±0.07 

9.16  ±0.35 

16 

4.5 

35 

Stello  et  al.  (2009) 

1.49  ±0.04 

42  ±  1 

19 

161797 

5317  ±  16 

5454  ±  35 

1.71  ±  0.01 

2.18  ±0.08 

28 

56.5  ±  0.07 

1200 

Bonanno  et  al.  (2008) 

0.87  ±  0.02 

951  ±  22 

21 

163588 

4451  ±  7 

4483  ±  25 

1 1.56  ±  0.06 

11.83  ±0.43 

2 

4 

36 

Stello  et  al.  (2009) 

1.36  ±0.03 

35  ±  1 

2 

168723 

4875  ±  7 

4923  ±  63 

5.92  ±  0.02 

7.43  ±  0.27 

26 

7.7 

80 

Barban  et  al.  (2004) 

0.68  ±  0.02 

64  ±2 

20 

181907 

5227  ±  199 

5637  ±  228 

12.10  ±  0.50 

12.01  ±0.52 

1 

3.47  ±0.12 

27b 

Carrier  et  al.  (2010) 

1.17  ±  0.17 

26  ±4 

5 

188512 

4992  ±  1 1 

5082  ±  69 

2.98  ±  0.01 

3.02  ±0.11 

1 

29.56  ±0.10 

472  ±  72 

Corsaro  et  al.  (2012) 

1.28  ±0.02 

470  ±9 

0.4 

Notes.  T\  is  the  interferometrically  calculated  effective  temperature  from  Table  2;  Tut  is  the  effective  temperature  from  the  literature  as  listed  in  Table  2.  Ri  is 
the  interferometrically  measured  radius;  Ra  is  the  radius  calculated  using  vmax  and  Av  from  asteroseismic  observations  as  well  as  7ia;  vmax  and  Av  are  from  the 
references  listed;  M  is  the  mass  calculated  using  Av  and  R\  \  and  vmax,caic  is  the  vmax  calculated  using  M  and  7) . 
a  No  vmax  was  listed,  so  it  was  calculated  using  Tcu  from  Martic  et  al.  (2001)  and  M  and  R  from  Kallinger  et  al.  (2009). 
b  No  vmax  was  listed,  so  it  was  calculated  using  Teff,  M,  and  R  from  Carrier  et  al.  (2010). 


HD  121370  and  HD  153210  display.  We  considered  whether 
or  not  this  would  have  an  effect  on  the  scatter  in  Figure  4  but 
the  stars  with  sparsely  sampled  curves  do  not  correspond  to  the 
outliers,  so  that  is  not  the  issue.  In  general,  we  do  not  observe 
any  systematic  trend  as  a  function  of  size.  The  residual  scatter 
in  the  giant  stars  is  comparable  to  what  Huber  et  al.  (2012b) 
found,  and  shows  the  relationships  between  observed  Av  and 
vmax  and  stellar  radii  are  not  as  precise  for  evolved  giant  stars  as 
they  are  for  dwarf  stars. 

The  scaling  relation  for  vmax  is  considered  to  be  less  robust 
than  the  relation  for  Av  (Huber  et  al.  2012b),  so  we  wanted  to 
test  it.  We  combined  Equation  (3)  with  our  R\  to  calculate  stel¬ 
lar  masses,  and  then  combined  the  masses  with  our  new  Te ff  to 
calculate  vmax  values  and  compare  them  to  the  measured  values. 
Table  4  lists  and  Figure  5  shows  the  results.  The  largest  outliers 
are  again  HD  153210,  HD  161797,  and  HD  168723  due  to  the 
discrepancies  in  calculated  radii  described  above.  HD  150680 
also  shows  a  16%  difference  between  the  observed  and  calcu¬ 
lated  vmax  value.  Our  angular  diameter  for  this  star  matches  those 
measured  by  Nordgren  et  al.  (2001)  and  Mozurkewich  et  al. 
(2003)  to  within  3a,  and  it  only  differs  from  the  SED  estimate 
by  3%.  We  note  that  the  mass  listed  in  Kallinger  et  al.  (2009)  is 
1.19  Mq  while  the  mass  determined  by  Martic  et  al.  (2001)  is 


1.3  to  1.5  Mq.  The  latter  agrees  with  the  mass  determined  using 
our  interferometric  radius  measurement:  1.33  ±  0.04  Mq.  In 
general  we  observe  good  agreement  between  the  observed  and 
calculated  vmax  within  the  uncertainty  bars  with  no  systematics 
with  respect  to  stellar  size  or  evolutionary  status. 

5.  SUMMARY 

We  measured  the  angular  diameters  of  10  stars  using  the 
NPOI.  The  combination  of  these  observations  with  other  infor¬ 
mation  from  the  literature  allowed  us  to  calculate  the  stars’  R , 
7’cn  ,  /"boi.,  and  L.  We  compared  our  interferometric  radius  R, 
values  to  those  determined  from  asteroseismic  scaling  relations 
and  found  good  agreement  between  the  two,  particularly  for  the 
less  evolved  stars.  Then  we  also  used  Av  from  the  literature  and 
our  /(,  to  calculate  the  stars’  masses  and  vmax  to  put  that  scaling 
relation  to  the  test  as  well.  Again,  the  results  agreed  to  within  a 
few  a  in  general. 

The  relations  work  best  for  main-sequence  stars  and  have 
limited  precision  for  giant  stars.  Hopefully  future  observations 
with  planned  spacecraft  such  as  Gaia  (Perryman  2003)  and 
TESS  ( Transiting  Exoplanet  Survey  Satellite ;  Ricker  et  al.  2009) 
as  well  as  planned  upgrades  to  existing  interferometers  such 
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Figure  4.  Comparison  between  interferometrically  measured  radii  and  those 
determined  asteroseismologically  listed  in  Table  4.  The  square  represents  the 
dwarf  star,  the  triangles  represent  subgiant  stars,  and  the  circles  represent  giant 
stars.  In  the  large  bottom  panel,  the  targets  from  Huber  et  al.  (2012b)  are  added 
in  as  stars.  The  dashed  line  is  the  1 : 1  ratio.  The  small  bottom  panels  show  the 
residuals  to  the  fit  normalized  to  the  asteroseismic  radii. 


as  the  addition  of  longer  baselines  on  the  NPOI  (increased 
resolution),  new  hardware  (increased  magnitude  limit),  and  the 
eventual  addition  of  large  telescopes  to  the  array  (see,  e.g., 
Armstrong  et  al.  2013)  will  lead  to  significant  improvements 
when  combining  data  from  two  techniques. 

The  Navy  Precision  Optical  Interferometer  is  a  joint  project  of 
the  Naval  Research  Laboratory  and  the  U.S.  Naval  Observatory, 
in  cooperation  with  Lowell  Observatory,  and  is  funded  by  the 
Office  of  Naval  Research  and  the  Oceanographer  of  the  Navy. 


Figure  5.  Comparison  between  the  calculated  and  observed  vmax  values  from 
Table  4.  The  symbols  are  the  same  as  in  Figure  4. 


This  research  has  made  use  of  the  SIMBAD  database,  operated 
at  CDS,  Strasbourg,  France.  This  publication  makes  use  of 
data  products  from  the  Two  Micron  All  Sky  Survey,  which 
is  a  joint  project  of  the  University  of  Massachusetts  and  the 
Infrared  Processing  and  Analysis  Center/California  Institute 
of  Technology,  funded  by  the  National  Aeronautics  and  Space 
Administration  and  the  National  Science  Foundation. 
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